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Abstract. Current status of dynamical modeling of relativistic heavy ion collisions 
and hydrodynamic description of the quark gluon plasma is reported. We find the 
hadronic rescattering effect plays an important role in interpretation of mass splitting 
pattern in the differential elliptic flow data observed at RHIC. To demonstrate this, 
we predict the elliptic flow parameter for <j) mesons to directly observe the flow just 
after hadronisation. We also discuss recent applications of outputs from hydrodynamic 
calculations to J/ip suppression, thermal photon radiation and heavy quark diffusion. 

1. Introduction 

r— i ; 

Understanding of the thermodynamic and transport properties of the quark gluon 
plasma (QGP) is one of the main topics in the physics of relativistic heavy ion collisions 
PQ. The QGP created in these reactions is by no means a static matter and exists as 
a transient state whose lifetime is of the order of ~ 1-10 fm/c. It is non-trivial even 
whether local equilibrium is really reached or not, which should be also clarified through 
the observation. Hence, one first has to justify the local equilibration by comparing the 
model calculation with experimental data. This is the starting point of further studies 
on the QGP under equilibrium. In order to understand thermodynamical (bulk) and 
transport properties of the QGP in relativistic heavy ion collisions, one then has to 
develop a dynamical analysis tool which links the observables with the equation of state 
(EoS) and transport coefficients of the QGP. It is inevitable to apply each relevant model 
to each stage of the collisions and unify these models consistently to describe the whole 
stage of the reactions. 

Along the line of thought, we have constructed a dynamical and unified model 
based on a fully three-dimensional hydrodynamic description of the QGP j2l [3] . One 
of the major discoveries at relativistic heavy ion collider (RHIC) is that elliptic flow 
parameters, which are intimately related with the degree of thermalisation, are well 
described within a framework of relativistic hydrodynamics [4] . An assumption of local 
thermalisation has to be made in hydrodynamics. However, local thermalisation can 
be expected only in the intermediate stage. Thus, one needs to model dynamics of 
reactions appropriately before and after the hydrodynamic stage to draw informations 



Hydrodynamic analysis of heavy ion collisions at RHIC 



2 



of the QGP from experimental data. Once we obtain the hydrodynamic evolution of 
the QGP, information of local thermodynamic variables is available to analyse other 
observables such as J/ip suppression, heavy quark diffusion and thermal radiation. In 
this way, hydrodynamic models enable one to analyse a vast body of data in relativistic 
heavy ion collisions. 

2. Dynamical modeling of heavy ion collisions 

Hydrodynamics is a general framework to describe the space-time evolution of locally 
thermalised matter for a given EoS. This framework has been applied to the intermediate 
stage in heavy ion collisions. We neglect the effects of dissipation and concentrate 
on discussion about ideal hydrodynamic models. The main ingredient in ideal 
hydrodynamic models is the EoS of hot and dense matter governed by QCD. In 
addition, one also needs to assign initial conditions to the hydrodynamic equations. 
Hydrodynamics can be applied to a system as long as its local thermalisation is 
maintained. However, in the final stage the particles are freely streaming toward the 
detectors and their mean free path is almost infinite. This is obviously beyond the 
applicability of hydrodynamics. Hence we also need a description to decouple the 
particles from the rest of the system. The hydrodynamic modeling of relativistic 
heavy ion collisions needs an EoS, initial conditions and a decoupling prescription. 
Hydrodynamic modeling has been sophisticated for these years and tested against RHIC 
data HI. 



The EoS is in principle calculated from lattice QCD simulations. Instead of using 
the result from lattice QCD simulations, we employ a phenomenological model in this 
study. At high energy density, the EoS can be described by a bag model 



Here P is pressure and e is energy density. The bag constant B is tuned to match 
pressure of the QGP phase to that of a hadron resonance gas at critical temperature 
T c : Pqgp(T c ) = Phadron(^c) ■ A hadron gas in heavy ion collisions is not in chemical 
equilibrium below the chemical freezeout temperature T ch which is close to T c [5], 
so the hadron phase may not be chemically equilibrated. A chemically frozen (or 
only partially equilibrated) hadron resonance gas can be described by introducing the 
chemical potential for each hadron [61 [7]. 

For a decoupling prescription, the Cooper-Frye formula [8] is almost a unique 
choice to convert the hydrodynamic picture to the particle picture. Contribution from 
resonance decays should be taken into account by applying decay kinematics to the 
outcome of the Cooper-Frye formula. The decoupling temperature T dec is fixed through 
simultaneous fitting of pt spectra for various hadrons in the low pt region. 

The prescription to calculate the momentum distribution as above is sometimes 
called the sudden freezeout model since the mean free path of the particles changes from 
zero (ideal fluid) to infinity (free streaming) within a thin layer. Although this model 
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is too simple, it has been used in hydrodynamic calculations for a long time. Recently 
one utilises hadronic cascade models to describe the gradual freezeout [21 El El [101 [11] 
in a more realistic way. As will be shown, the hadronic afterburner is mandatory in 
understanding elliptic flow data. Phase space distributions for hadrons are initialized 
below T c by using the Cooper-Frye formula. We employ a hadronic cascade model JAM 
[T2] to describe the space-time evolution of the interacting hadron gas. This kind of 
hybrid approaches in which the QGP fluids are followed by hadronic cascade models 
automatically describes both the chemical and thermal freezeout and is much more 
realistic especially for the late stage than the naive Cooper-Frye prescription. 

Initial conditions in hydrodynamic simulations are so chosen as to reproduce the 
centrality and rapidity dependences of multiplicity dN c ^/drj. Initial conditions here are 
energy density distribution e(x,y,r] s ) and flow velocity y, r) s ) at the initial time 
To. Again baryon density is neglected since the net baryon density is quite small at 
midrapidity at RHIC. At mid-(space-time)rapidity r/ s = 0, one can parametrise the 
initial entropy density based on the Glauber model 



The soft/hard fraction a is adjusted to reproduce the measured centrality dependence 
[T3] of the charged hadron multiplicity at midrapidity. By using the EoS, one can 
calculate the initial energy density distribution together with pressure and temperature 
distributions from Eq. (T5]). For space-time rapidity dependence of initial conditions 
and a novel initial condition based on the colour glass condensate (CGC) picture, see 
Refs. [21 [H]. The fitting of multiplicity is the starting point of further analysis based 
on hydrodynamic simulations. 

In the hydrodynamic models, various combinations of initial conditions, EoS and 
decoupling prescriptions are available to analyse the experimental data. Of course, final 
results largely depend on modeling of each stage. So it is quite important to constrain 
each model and its inherent parameters through systematic analysis of the data toward 
a comprehensive understanding of the QGP. 

3. Elliptic flow 

Anisotropic transverse flow can be quantified by Fourier components of azimuthal 
distribution, v n = (cosny?), where the average is taken over by weighting azimuthal 
distribution. The second harmonics t>2 is called elliptic flow parameter and is intimately 
related with the EoS and the degree of thermalisation |15j . 

3.1. Charged hadrons 

The pseudorapidity dependence of v 2 observed by PHOBOS [16] has a triangular shape: 
V2 is maximum at midrapidity and decreases as moving away from midrapidity. On the 
other hand, in the ideal hydrodynamic calculations with T dec = 100 MeV, t>2 depends 
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Figure 1. (Left) Pseudorapidity dependence of V2 for charged hadrons. PHOBOS data 
[16] are compared to the ideal hydrodynamic model with T dcc = 100 and 169 MeV 
and the hybrid model [2]. (Right) V2 for charged hadrons as a function of centrality 
in | ?y |< 1. PHOBOS data [H] are compared to the ideal hydrodynamic model with 
T doc = 100 MeV and the hybrid model 0. 

mildly on pseudorapidity in | rj |< 3. Consequently, the hydrodynamic model generates 
as large v 2 as the data only near midrapidity and overshoots the data significantly at 
forward and backward rapidities j?l [T7j as shown in Fig. [I] (left). If we replace the 
hadron fluid with a hadron gas utilising a hadronic cascade, V2 is generated little in 
the forward and backward region during hadronic rescattering stage [2j. In this hybrid 
model the hadrons have a finite mean free path and exhibit an effective shear viscosity 
in the hadron phase. So the dissipative hadronic "corona" effect [H] turns out to be 
important in understanding the charged hadron t> 2 data. Figure [1] (right) shows the 
centrality dependence of t>2 at midrapidity (| 7] |< 1). The solid line is the result from 
an ideal hydrodynamic calculation, while the dotted line from the hybrid model. It is 
clear that for peripheral collisions where the multiplicity is small the hadronic viscosity 
plays an important role in description of the t>2 data. One may notice that the result 
from the hybrid model systematically and slightly smaller than the data. However, 
there could exist an effect of initial eccentricity fluctuations which is absent in the 
present hydrodynamic calculations and could enhance t> 2 by ~ 90% in 0-5% and by 
~ 20% in 50-60% centralities [19]. The deviation between the results and the data may 
be interpreted quantitatively by this eccentricity fluctuation effect. 

3.2. Identified hadrons 

The hybrid model also reproduces a mass ordering pattern of V2 for identified hadrons 
as a function of near midrapidity [20J (Fig. [2j left) and in forward rapidity region 
[2~T] . We consider semi-central collisions (20-30% centrality), choosing impact parameter 
b = 7.2 fm. If one would look at the result just after the hadronisation, the difference 
between pions and protons would be quite small as shown in Fig. [2] (right). As 
demonstrated in this figure, at low the pion curve moves up while the proton curve 
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Figure 2. (Left) V2(pt) for identified hadrons. STAR data [20] are compared to the 
hybrid model [3j. (Right) V2(pt) with and without hadronic rescattering [3]. 



moves down through rescatterings. So the hadronic rescatterings are attributed to 
the observed splitting patterns of ^(pr)- It should be noted that the same hadronic 
rescattering effect is also seen in the recent study based on the different hadronic 
cascade model [22J. One can conclude that the large magnitude of the integrated Vi 
and the strong mass ordering of the differential v^ipr) observed at RHIC result from a 
subtle interplay between perfect fluid dynamics of the early QGP stage and dissipative 
dynamics of the late hadronic stage. The large magnitude of V2 is due to the large overall 
momentum anisotropy, generated predominantly in the early QGP stage. Whereas the 
strong mass splitting behavior at low pt reflects the redistribution of this momentum 
anisotropy among the different hadron species during the hadronic rescattering phase. 

Hadronic rescattering turns out to cause the mass ordering of V2(pr)- Then, what 
happens to the (hidden) strangeness sector in which rescattering with other hadrons 
does not happen so frequently [23]? Figure [3] (left) shows v 2 (pt) for ir, p and 0. The 
V2(pt) curves for pions and protons separate as discussed before, whereas V2(pt) for the <fi 
meson remains almost unchanged during the hadronic stage. As a result of rescattering 
the proton elliptic flow ends up being smaller than that of the <p meson, v^ipr) < vt(Pr) 
for <Pt < 1-2 GeV/c, even though > m p . Hadronic dissipative effects are seen to 
be particle specific, depending on their scattering cross sections which couple them to 
the medium. The large difference of cross section between protons and 4> mesons in the 
hadronic rescattering phase leads to a violation of the hydrodynamic mass ordering at 
low pt in the final state. Note that clear mass ordering is seen in Fig. [3] (right) if (j) 
mesons were also assumed to fully participate in hydrodynamic flow. 

4. Application of hydrodynamic results 

So far, we have seen that the hybrid model, in which the dynamical evolution of the 
locally thermalised QGP and the dissipative hadronic gas are incorporated, works well 
in describing the elliptic flow data. One can utilise the space-time evolution of the QGP 
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Figure 3. V2(pr) for pions (dotted), protons (dashed) and 4> mesons (solid) in Au+Au 
collisions at b = 7.2 fm from the hybrid model (left) and from the ideal hydrodynamic 
model with T dcc = 100 MeV (right). 



obtained above to study other observables. 

The charmonium is one of the major diagnostic tools of the QGP in relativistic 
heavy ion collisions [21] . The melting temperature of J/ip, if)' and \c is obtained by 
modeling the free streaming of charmonia on the QGP background [251 126] . We first 
generate the charmonia in the transverse plane, distribute them according to the binary 
collision density, simulate the transverse motion of charmonium and then dissociate them 
when the temperature at the position of a J/if> (if)' and Xc) is larger than the melting 
temperature Tj/^ (T x ). It is found [251 12S] that the suppression pattern as a function 
of centrality reflects the region of the temperature above Tj/^ which is obtained in 
hydrodynamic calculations. The best fit to the suppression factor observed by PHENIX 
[27] is achieved by choosing Tj/^ = 2T C and T x = 1.34T C . 

Thermal photon is one of the direct and penetrating probes to study the QGP 
[28J . To calculate the transverse momentum spectra of thermal photons, one convolutes 
the reaction rate, which is the yield per unit space-time volume, over the volume 
of thermalised matter in heavy ion collisions which is obtained from hydrodynamic 
simulations. By using the hydrodynamic results obtained in the previous section, we 
calculate the thermal photon spectra [29j [30] as well as components of hard photons, 
jet-photon conversion and fragmentation. We find that thermal photons originated from 
the QGP phase are dominant around pr ~ 3 GeV/c [29, 30J. 

From hydrodynamic simulations in the previous section, the maximum temperature 
is found to be at most ~400 MeV in central collisions which is much smaller than the 
mass of heavy quarks. So it is safe to assume that these heavy quarks are not thermally 
excited and are produced only in the first hard collisions. We treat these heavy particles 
as "impurities" in the medium in which light degrees of freedom such as u, d and s 
quarks and gluons are equilibrated. We model the transverse motion of heavy quarks 
as the Brownian motion and calculate the nuclear modification factor for electrons and 
positrons from semileptonic decays of heavy mesons and elliptic flow parameter for heavy 
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Figure 4. (Left) Nuclear modification factor for electrons and positrons from 
semileptonic decays in semi-central collisions for 7 = 0.3, 1.0 and 3.0. Experimental 
data are obtained by PHENIX 32J. Upper and lower bands represent when diffusion 
of heavy quarks stops in the mixed phase /qgp = 1 an d 0, respectively, where /qgp 
is the fraction of the QGP phase in the mixed phase. The middle points in the band 
correspond to /qgp = 0.5. (Right) Sensitivity of v 2 (pt) for bottom quarks to the 
parameter in the drag force. Results arc obtained with /qgp = 0.5. 



quarks. We solve the following relativistic Langevin equations 

M(t) = ^-rAt, Ap(t) = -rpAt + g{t) (3) 

on the hydrodynamic background |31j . Here M, x and p are mass, position and 
momentum of a sample heavy quark. £ is the Gaussian white noise. We parametrize 
drag force as T = jj?. The above equations are solved in the local rest frame of a 
fluid element which is moving with in the center of mass system. 7 in the drag force 
is the only adjustable parameter in this modeling of heavy quark diffusion. From a 
comparison of nuclear modification factor with the data in Fig. H] (left), we find 7 ~ 1-3 
gives a reasonable description of the data observed by PHENIX [32]. If we look at V2(pt) 
for heavy quarks just before hadronisation, the results with 7 ~ 1-3 are considerably 
smaller than the result with 7 = 30.0 whose relaxation time r = 1/T is much smaller 
than the life time of the QGP. This indicates that the heavy quarks do not reach the 
complete thermalisation limit. 



5. Summary 

We discussed the current status of hydrodynamic description of relativistic heavy ion 
collisions. Integrated and differential elliptic flow parameters in low p? regions are well 
reproduced by the hybrid approach in which ideal hydrodynamic description of the 
QGP is followed by microscopic description of hadron gas. We found the observed mass 
ordering pattern results from the hadronic rescatterings. As a consequence, violation 
of mass ordering is predicted for mesons that do not couple to the hadronic medium. 
We also discussed the recent application of hydrodynamic results to J / ip suppression, 
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thermal photon radiation and heavy quark diffusion. 
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